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Evolution of a particle in an inverse square potential is studied. We 
derive an equation of motion for {r'^'^ and solve it exactly. It gives 
us a possibility to identify the conditions under which a falling of a 
quantum particle into an attractive centre is possible. We get the time 
of falling of a particle from an initial state into the centre. An example 
of a quasi-stationary state which evolves with being constant 
in time is given. We demonstrate the existence of quantum limit of 
falling, namely, a particle does not fall into the attractive centre, when 
coupling constant is smaller then some critical value. Our results are 
compared with experimental measurements of neutral atoms falling in 
the electric field of a charged wire. Moreover, we propose modifications 
of the experiment, which allow to observe quantum limit of falling. 
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1 Introduction 

The attractive inverse square potential is studied in literature from different 
points of view. Such interest is motivated by peculiarities of this potential. 
There are no stationary energy levels for the attractive 1/r^ potential, an 
energy of a particle in this potential is not bounded from below [DEiia]. 
Formally, one can obtain solutions of the Schrodinger equation 

ih ( 1 ) 

as Bessel functions Z,,{kr), where u = ^1/4 + l^/h? — 2m'y/h?, k = ^J2mE/h 
[1]. The corresponding functions fjE are integrable at the inhnity, but in the 
vicinity of the attractive centre they have asymptotic fjE ~ + Br^-, 

where s± = —1/2 ± ^Jl/A — 2m'y / h ?, and are nonanalytical for a strong 
coupling regime 7 > h^/8m. The nonanalytical behaviour of 'ipE can be as¬ 
sociated to the classical falling into the centre |3]. In addition, Hamiltonian 
(IT!) is not hermitian [2]. 
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Despite the fact that —1/r^ potential is quite exotic, it appears in such 
physical problems as the Ehmov effect [1], a neutral atom in the electric held 
of a thin charged wire EE], an atom with magnetic moment in the magnetic 
held of a long solenoid [7j, the matter near horizon of a black hole [8], an 
electron near a dipolar molecule [9]. 

Bound states can be restored in an inverse square potential by building 
self-adjoint extensions of ([I]) [in],[I2]- Also bound states for an inverse square 
potential may appear within the renormalization scheme ini.ini- It is in¬ 
teresting that an inverse square potential is regularized in a natural way in 
a space with a minimal length na and noncommutative space IH- In all 
these methods an additional parameter with the dimension of length appears, 
which is not a parameter of 1/r^ potential. 

The solution of the classical problem is known for a long time, from the 
early 18th century [T5|. Trajectory of a particle moving in an inverse square 
potential is one of the three possible Cotes spirals. A particle can either fall 
into the attractive centre or escape from it. The situation depends on the 
value of a coupling constant, z-component of the angular momentum, and 
initial position. 

The aim of this paper is to study the evolution of a quantum particle in 
an inverse square potential. On the contrary to the classical case, evolution 
of a quantum particle in an inverse square potential has not been discussed 
in literature. It is an interesting question how does a wave function evolve 
under Hamiltonian ([T]), how does a falling into the centre look for a quantum 
particle. 

Our paper is organized as follows. In Section 2, we consider the evolution 
of (r^) in an inverse square potential. Our theoretical results are compared 
with experimental measurements of falling of Lithium atoms into the charged 
wire [B] in Section 3. In Section 4, we come up with the description of 
modihcations of this experiment, which allows to observe the quantum falling 
limit experimentally. We summarize our results in the concluding section. 

2 A time of falling 

Calculating an evolution of a wave function under a Hamiltonian with dis¬ 
crete spectrum can be easily done by expanding the wave functions in eigen¬ 
states of the Hamiltonian. Doing so, one can act trivially on this expansion 
by the evolution operator and obtain how the wave function changes over 
time. But in the case of an inverse square potential, where no stationary 
levels exist, this scheme cannot be performed. 

Another way of calculating the evolution is to expand the evolution opera- 


2 


tor into the Taylor series by time. For small times this series might converge. 
But for inverse square potential in this scheme additional problem appears 
to be related to the singularity of the Hamiltonian. To obtain the evolution 
of the wave function, one should analyse all terms in Taylor expansion. If 
the wave function behaves as near the origin, then the n-th order term of 
expansion behaves as and is singular for n > s/2. 

The evolution of operator averages can be obtained from Heisenberg equa¬ 
tions. Let us consider the evolution of (r^). This average reflects a spatial 
distribution of the wave function, namely a smaller (r^) corresponds to a 
more compact localization of the wave function, and vice versa. A particle 
collapses into a point when (r^) = 0. The hrst derivative of (r^) with respect 
to time reads 


/L 

dt 
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The second derivative gives 
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If H is a power function in r, one obtain a chain of equations, which can be 
closed only for three potentials, namely H = 0, V = kr"^, and V = 

In the case oiV = —y/r^ we hnd 






( 4 ) 


Since the energy is a constant of motion, the solution of (jl]) can be im¬ 
mediately found in the following form 


+ tE+irkt + ( 5 ) 

where {...)o denotes averaging at the initial moment of time f = 0. 

Equation of evolution for (r^) (jS]) coincides with the corresponding clas¬ 
sical equation where operator averages are replaced by their classical values. 
Classical equation can be rewritten in the following manner 

= rl + 2rorot H- f. (6) 

m 

The fate of a particle depends on the values of (r^)g, {rp + pr)^, and (H). 
In the case of (H) < 0 a particle obviously falls into the centre. For (H) > 0 
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a particle either escapes from the attractive centre or falls into the centre. 
The falling is possible if {rp + pr)^^ < 0 and (r^)^ < (rp + pr)^/8m |(i7)|. 
Ineqnality (rp + pr)^ < 0 corresponds to the classical case of tq < 0, i.e. 
radial projection of the initial velocity is directed to the force centre. An 
analogous situation is in the case of {H) = 0. It is interesting that for 
{H) = 0 and {rp + pr)^ = 0 the particle moves with a constant average 
(r^). The time dependence of (r^) is plotted in the Figure 1. Falling time 
can be calculated from ([5]) by putting (r^) = 0. 
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FIG. 1. Evolution of (r^) for different initial conditions. The complete line 
(dashed and solid) represents the solution of the equation of motion for (r^), 
but a real motion of a particle is pictured by the solid part. The time scale fac¬ 
tor to is defined in the following way. If {H) ^ 0, then Iq = ^Jm /2 |(Fr)| 
and e = {rp pr)^ / ^J2m \{H) \ (r^)^. In the case of {H) = 0 we can put 
to = m (r^)Q / |(rp + pr)^]. If both {H) and (rp -|- pr)^ equal zero, particle 
evolve with a constant (r^), and in this case the unit of time can be chosen 
arbitrarily. 
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In the next section we will compare our theoretical results with measure¬ 
ments of neutral atoms motion in the electric held of a charged wire [B] (from 
here and through the text we will call it as DUS experiment due to names 
of authors). In this experiment particles move in two-dimensional potential. 
So, the 2D case will be considered in detail. Due to the operator identity 
rp + pr = r-^ + -^r, in 2D space when the radial part of the wave function 
is real we have 


{rp + pr)^ = 0. (7) 

Taking into account condition ([7]), we can obtain time of falling as follows 


tf = 



( 8 ) 


If dZ]) holds, the particle falls into the centre only if its energy is negative. 
If the wave function is proportional to near the origin with s < 1 / 2 . (H) 
formally becomes inhnite, because corresponding integrand is singular. In 
this case, the time of falling equals zero. But such states require an inhnite 
energy for creation, therefore they are unphysical. For positive energies tf is 
imaginary, and particle cannot fall into the centre. 

It is interesting that a particle with zero energy will stay inhnitely long 
on the same distance from the centre (r^) = (r^)g; this fact is easily seen 
from (|5]). In this case, the time of falling formally is inhnite. An example of 
such a quasi-stationary state can be found in the following form 


i’s 


2/3*+i 

r(s + i)’ 


' exp(' 



(9) 


where r(s) is the Euler’s gamma function. The energy is equal to zero for 
s = 2 m 7 / 3 h^. This condition provides the quasi-stationary motion. 

Classical expression for a falling time can be obtained from the quantum 
one ([H]) by replacing operator averages by their classical values. The principal 
diherence between classical and quantum evolution is in a set of allowed val¬ 
ues of classical quantities and corresponding quantum mechanical averages. 
An interesting quantum ehect appears when 7 is small. The classical particle 
can fall into the centre for arbitrarily small 7. If one put r = 0 and L = 0 
then E becomes negative, the particle will fall into the centre. In quantum 
mechanics situation is quite diherent, there is a critical value 7 c, so that 
for 7 < 7c particle will not fall into the centre. This quantum limit exists 
due to the fact that for small 7 an energy of a particle cannot be negative. 
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To obtain this critical value 7 c let us analyse the average of Hamiltonian 
in detail. For the sake of simplicity, we take the wave function in the form 
i) = 7^(r)e*'"‘^/\/^, where radial part TZ{r) is an arbitrary function, the 
angular part is an eigenfunction of orbital momentum, which corresponds to 
the eigenvalue Iz- In this case one can rewrite (H) in the following form 




2m 2m 





where pr = ■\dr\/r. Using the condition {H) < 0, we get 


7 > 7c 


{pD , 

2m(l/r2) 2m 



( 10 ) 


( 11 ) 


From flTT]) it follows that min 7 c = Fc = ^. It is remarkable that Fc 
coincides with the limit of a strong coupling regime. For a given wave function 
the critical value 7 c can be calculated from flTT]) . For example, for the wave 
function 'ipg dH]) (integrand is nonsingular when s > 1/2) we obtain 


s + l f P? 7 \ 

(r2)p \2m s) ’ 

and for the critical value 7 c we get 



sh^ 


4sFc. 


( 12 ) 


(13) 


3 Neutral atoms in the field of a thin charged 
wire 

Falling into an inverse square potential was experimentally realized in DUS 
experiments with ultracold Lithium atoms in the held of a thin charged wire 
[6]. In this experiment neutral Lithium atoms with polarizability au were 
loaded in a cylindrical chamber with length of Z = 0.1 m. The thin wire with 
a radius of ri =0.7 pm was placed coaxially with the chamber. Radius of the 
external cylinder is not pointed out in papers, but it can be obtained using 
the relation between the charge on the wire and the voltage on it (external 
cylinder is grounded): a charge of g = 640 ^ corresponds to the voltage 
of U = 100 V. Using the expression for the capacitance of a cylindrical 
capacitor, we obtain the chamber radius as r 2 = 4.2 mm. In this experiment 
the authors achieved the pressure of p = 6 ■ 10“^° Torr. This pressure is 
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sufficiently small, thereby we can assume that Lithium atoms do not collide 
with each other during the motion. 

Neutral atoms are polarized by the electric held of a wire and obtain 
an electric dipole moment SoaE. The Hamiltonian of the induced electrical 
dipole in an electric held E reads 

= p 1 _ £22L! = ih _ g, (14) 

2m 2 2m 

where 7 = ag^/87r^eo- 

Now we want to estimate the time of falling of atoms in this chamber using 
expression ([ 8 ]). Let us assume that the wave function of the Lithium atoms 
is equal to ([9]) with s = 1/2. Putting (r^) = and using (jS]) and flT^ . we 
obtain the time of falling of about ~ 7 s. This result is in correspondence 
with the estimation made in [B], where the lifetime of atoms of about 10 s 
was measured. 

An exact wave function of Lithium atoms is not known. So, to estimate 
a credible value of critical charge Qc needed for a falling of quantum particle 
we have to use expression for T,.. Finally, we obtain 


Qc 


h^Tl'^er 


ma 


( 15 ) 


In the classical limit {h —)■ 0) the critical charge vanishes Qc = 0, as it should 
be. 

The critical charge of the wire in DUS experiment, for which Lithium 
atoms {aLi = 24.3 A ) do not fall into the centre, is about a g = 1.8 —. 
This charge corresponds to a voltage of Uc = 0.29 V. Because of a noise 
produced by optical trap in the chamber such value of voltage cannot be 
measured in this experiment. 


4 Proposed experiment for observing quan¬ 
tum limit 

The aim of this section is to give a proposal how to increase the critical 
value Uc and make it observable. The hrst proposal is to use lighter atoms 
with smaller value of polarizability, as it follows from expression for Qc (na. 
For Hydrogen atoms {an = 0.667 A^) the critical charge and voltage for the 
chamber in the DUS experiment are Qc = 30 Uc = 4.6 V respectively. 
So, for Hydrogen atoms values of quantum limits are 16 times bigger, than 
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for Lithium atoms. Almost the same values are obtained for ^He atoms: 
g = 31 2^, f/c = 4.8 V. 

Another idea is to change a geometry of the chamber. This proposal is 
motivated by the relation between the charge and the voltage of capacitor 


U = 


a 

C 


^In^, 

27r£o n 


(16) 


where C is the capacitance. So, using a capacitor with bigger r 2 and smaller 
ri allows us to work with bigger voltage. For another capacitor with = 
X 2 r 2 , r[ = Airi the value of the critical voltage reads 


U' = U. 


In^ 

Airi 


In Aa 


(17) 


For DUS experiment we have In — 8.67, so the critical voltage is f/' ~ 

1.26Uc for A 2 /A 1 = 10 (r 2 = 4.2 cm, r 2 = 0.7 pm) and U' 1.5374 for 
A 2 /A 1 = 100 (r 2 = 42 cm, r 2 = 0.7 pm). Increasing the size of chamber by 
10 times ensures the critical voltage of about 6.1 V for ^He atoms and 5.8 V 
for H atoms. It is almost 21 times bigger voltage than the critical voltage 
in the original DUS experiment. We hope that such improvements give a 
possibility to observe the quantum limit of falling experimentally. 


5 Conclusions 

We have shown that an average (r^) evolves as a quadratic polynomial in time 
in an inverse square potential. This evolution is fully dehned by averages of 
(r^)Q, {rp + pr)^, and {H) in the initial state and is qualitatively different 
for different initial conditions (Fig. 1). We have found necessary conditions 
needed for a falling of a quantum particle into the attractive centre. The case 
of {rp + pr)^ = 0 in two-dimensional space has been considered in detail. 
For such states the time of falling is given by expression ([8]). We have shown 
that there are quasi-stationary states, which evolve with constant (r^). Also 
an example of such a state has been given ([9]). By analysing the average (77), 
we have established a quantum limit of falling, namely a particle cannot fall 
into the attractive centre if the coupling constant 7 is smaller than some 
critical value 7 c. The critical value 7 c is dehned by the expression (ITT]) . 

We have compared the obtained results with the experimental measure¬ 
ments of the motion of neutral Lithium atoms in the electric held of a charged 
wire |6]. Calculated theoretically time of falling coincides with the experi¬ 
mentally measured value. Also, we have calculated the critical charge of the 





wire and the voltage between the wire and the wall of a chamber, which 
allows qnantnm particle to fall into the centre. The critical voltage is abont 
0.3 V. 

Unfortunately, critical parameters of quantum falling of Lithium atoms 
in DUS experiment are so small that they could not be measured there. 
We have given some proposals how to improve the experiment to observe 
the quantum falling limit. Namely, we propose to use lighter atoms with 
smaller polarizability (Hydrogen or Helium atoms) and to change the size of 
a chamber. This improvements make a voltage about 6.1 V, that is almost 21 
times bigger than analogous values in DUS experiment. We hope that these 
improvements allow to observe the quantum limit of falling in the experiment. 
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